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i-^ ■ Abstract 

> ' 

• I— I . We study analytically the phase diagram of the pure SU{N) lattice gauge 

K> ' theory at finite temperature, and we attempt to estimate the critical decon- 

C^ ■ finement temperature. We apply large A^ techniques to the Wilson and to 

the Heat Kernel action, and we study the resulting models both in the strong 
coupling and in the weak coupling limits. Using the Heat Kernel action, 
we establish an interesting connection between the Douglas-Kazakov phase 
transition of two-dimensional QCD and the deconfining transition in d di- 
mensions. The analytic results obtained for the critical temperature compare 
well with Montecarlo simulations of the full theory in (2 -|- 1) and in (3 -|- 1) 
dimensions. 
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1 Introduction 

Finding a precise characterization of the finite temperature deconfining transition 
is a long standing problem of Lattice Gauge Theory (LGT). Since the seminal work 
of Svetitsky and Yaffe [jl|, it is understood that all the relevant properties of the 
transition can be encoded in a suitable effective action for the order parameter, the 
Polyakov loop, and a qualitatively satisfactory description of the phase diagram 
can be obtained in this way. Such an approach clearly involves an enormous sim- 
plification of the model, through a drastic reduction of its degrees of freedom; the 
non-trivial task is then the construction of an appropriate effective action. The 
simplest and most popular proposal for such an action is 0-0 

Seff{J) = Y: Re\jY: TT,{V{x))TT,iV\x + z))\ , (1) 

where Trf(V"(x)) is the Polyakov loop. 

This action has been widely studied in the literature, both at finite A^ [^ |^ and in 
the large N limit [^, Q , with strong coupling p, ^ and mean field approximations, 
and with Montecarlo simulations [@, §• In particular, in the large N limit it is 
exactly solvable, and leads to a peculiar phase diagram, with a first order point and 
a line of third order transitions of the Gross- Witten type [^ ^. In the literature, the 
action (|I]) has been mainly used as a tool to understand the general features of the 
phase diagram, without really hoping to obtain precise determinations of, say, the 
location of the critical deconfinement temperature T^. This is due to the fact that 
the approximations involved in going from the complete LGT to eq. (P are rather 
drastic, and the correspondence between the coupling J in eq. (|l]) and the coupling 
(3 of the complete LGT (and the consequent possibility of finding a meaningful 
result for the critical temperature in the continuum limit) is rather ambiguous. 

One of the goals of the present paper is to show that, by using new techniques 
which have become available in these last years for the large N analysis of gauge 
theories, one can go a bit further and extract from a suitable generalization of the 
action (|ID some reliable information on the deconfinement temperature. Even if 
the values obtained are still affected by large uncertainties, they essentially agree 
with those extracted from Montecarlo simulations. This is a significant progress, 
because the results are obtained analitycally and, what is more important, our 
approach opens the way to further improvement. In fact a better control on the 
approximations involved may be obtained within the framework of the hot Eguchi- 
Kawai model, and also the extension to the case of finite N may be studied. 

In this paper, we shall show that the action (|I|) can be considered as an ap- 
proximation of a much more accurate effective action which, despite its complexity, 
can still be precisely studied in the large A^ limit, in weak coupling as well as in 
strong coupling. Besides describing the phase diagram of the model, we will make 
an effort to extract as precisely as possible the value of the critical temperature 
in terms of the coupling of the original model, studying its scaling behaviour and 



making contact with the known results from Montecarlo simulations both in (2+1) 
and (3+1) dimensions. As we shall see in the last section, the two cases require 
very different analysis, and only in the (2+l)-dimensional case the approximations 
are really under control, the value that we obtain has the correct scaling behaviour, 
and the continuum limit can be correctly approached. In (3+1) dimensions things 
are significantly more complicated, and our results should be better taken as strong 
coupling estimates, valid for small values of the Nf (the lattice size in the time direc- 
tion); nevertheless, also in this case, some interesting information can be extracted, 
and a rather succesfuU comparison with known Montecarlo estimates can be made. 
To describe QCD on the lattice, we will use both the Wilson action and the 
action defined in terms of the Heat Kernel on the group manifold [^ , already em- 
ployed in this context by |]T0| . We will study the approximations that are needed to 



reduce these actions to the model (|I|), namely a Migdal-Kadanoff renormalization 
and the decoupling of the space-like plaquettes. We observe that the Heat Ker- 
nel action has definite advantages from the point of view of analytic calculations, 
and, upon decoupling the space-like plaquettes, allows an exact integration over the 
space-like degrees of freedom. The resulting model can be studied, in the large N 
limit, both in weak coupling (where it is directly equivalent to the Wilson action), 
and in strong coupling, where a non-trivial transformation is needed to make the 
connection with the Wilson action. 

In the process, we shall make contact with some interesting result recently ob- 
tained in the context of two-dimensional QCD (QCD2). In particular, we shall 
discuss the relation between the deconfinement transition in {d + 1) LGT and the 



so called Douglas-Kazakov transition [11] in QCD2. 



Our analysis starts with a review of the various approximations needed to obtain 
Seff- This will be done in the next section, where we shall also fix our notations 
in the general framework of finite temperature LGT. Section 3 will be devoted to 
a discussion of how the same approximations, applied to the Heat Kernel, link the 
effective action for the Polyakov loop to the two-dimensional QCD on a cylinder. 
We then review some of the features of QCD2 which are relevant in this case. 
Sections 4 and 5 will be devoted to a weak coupling and a strong coupling expansion 
respectively, while in the last section we shall make contact with the known results 
on the deconfinement temperature of LGT, and we shall make some concluding 
remarks. 



2 Finite Temperature LGT 

Let us consider a pure gauge theory with gauge group SU{N), defined on a rf + 1 
dimensional cubic lattice. In order to describe a finite temperature LGT, we have to 
impose periodic boundary conditions in one direction (which we shall call from now 
on "time-like" direction), while the boundary conditions in the other d direction 
(which we shall call "space-like") can be chosen freely. We take a lattice of A^j (A^^) 



spacings in the time (space) direction, and we work with the pure gauge theory, 
containing only gauge fields described by the link variables Un-i G SU{N), where 
n = {x, t) denotes the space-time position of the link and i its direction. It is useful 
to choose different bare couplings in the time and space directions. Let us call them 
/3t and (3s respectively. The Wilson action is then 

Sw = T. J^^Af^tT, Trf (t/„;o^) + /^s E Trf (?7„,,) I , (2) 

where Trf denotes the trace in the fundamental representation and Un-oi {Un;ij) are 
the time-like (space-like) plaquette variables, defined as usual by 

U^,, = Un,U^^,,Ul^.^Zli, . (3) 

In the following we shall call Ss (St) the space-like (time-like) part of Sw- Ps 
and Pt are related to the (bare) gauge coupling g and to the temperature T by the 
usual relations 

2N 



where a is the space-like lattice spacing, while jpj; is the time-like spacing. The two 
are related by the dimensionless ratio e = jpfa- We can solve the above equations 
in terms of e, as 

A = ^a^-^ , (5) 

re 

Ps=^-^a'-' . (6) 

In a finite temperature discretization it is possible to define gauge invariant 
observables which are topologically non-trivial, as a consequence of the periodic 
boundary conditions in the time directions. The simplest choice is the Polyakov 
loop, defined in terms of link variables as 

Nt 

P{x) = TT,V{x)=TT,ll{U,,t,o) ■ (7) 

t=i 

As it is well known, the finite temperature theory has a new global symmetry 
(unrelated to the gauge symmetry), with symmetry group the center C of the gauge 
group (in our case Z^). The Polyakov loop is a natural order parameter for this 
symmetry. 

In (i > 1, finite temperature gauge theories admit a deconfinement transition at 
T = Tci separating the high temperature, deconfined, phase (T > Tc) from the low 
temperature, confining domain (T < Tc). In the following we shall be interested in 
the phase diagram of the model as a function of T, and shall make some attempt 



to locate the critical point Tc. The high temperature regime is characterized by the 
breaking of the global symmetry with respect of the center of the group. In this 
phase the Polyakov loop has a non-zero expectation value, and it is an element of 
the center of the gauge group (see for instance |jl|). 

The Wilson action eq.(|^) is clearly too complex to be handled exactly, and some 
approximate effective action is needed. As we shall see, at least three different 
approximations are needed, in order to reach a solvable model. The first step is to 
identify the Polyakov loops as relevant dynamical variables. This can be done by 
dimensional reduction, as described in the next subsection (2.1). Second, in order 
to integrate out the spacelike degrees of freedom, one must decouple the spacelike 
plaquettes, and this will be described in subsection (2.2). The third and last step 
is to take the large A^ limit, which allows to solve exactly the model, through a 
large A^ factorization of the Polyakov loops (subsection 2.3); the large A^ limit is 
important even when it does not lead to exact solvability, because in that limit a 
translational invariant master field is expected, and the model can be reduced to 
an effective one-plaquette model. 

In the remaining part of this section we shall describe in some more detail 
these steps, and discuss their reliability. Since the approximations are well known, 
and have been already widely discussed in the literature, we shall mostly state the 
results, and refer the reader to some of the original papers. 

2.1 Migdal-KadanofF approximation 

By means of an approximate renormalization group transformation, within the 
framework of a Migdal-Kadanoff bond-moving scheme, it is possible to reduce the 
original {d + 1) dimensional Wilson action, eq. (|^) to a rf dimensional gauge the- 
ory coupled to the Polyakov lines, which play the role of a Higgs field p|, |^. The 
resulting action is 



exp(S'i 



n 




j=i 



X exp 



iNt 



drXriVix)U£.,V\x + t)Ul.,) 



X 



-i?e /3.$: Trf(f/,.,,, 



Here the sum is limited to the d spacelike directions, Ir{Pt) is the coefficient of 
representation r in the character expansion of the Wilson action, and the new 
spacelike coupling is related to the old one by 



Ps = Nt/3s 



(9) 



In the literature the timelike part of this action is then further approximated by 
truncating the sum over representations at the first term, so as to reconstruct once 
more a Wilson-type action 



^2 = E T^ ^e i A X: Tt,{V{x)U,.,V^{x + z)UI^,) + 4 ^ Trf(f/.-.,) I , (10) 

X [ i=l i<j J 

As we shall see, this truncation is not necessary if one works with the Heat 
Kernel action, to be discussed in the next section. In the present framework, the 
Migdal-Kadanoff approximation is expected to be a good approximation both at 
strong and at weak coupling (see for instance [0)) while it might be inadequate at 
intermediate couplings. In particular, it is easy to see that the timelike part of 5*1, 
eq. (^, could have been obtained also within the framework of a strong coupling 
expansion. 

The truncation of eq. (^ to eq. ( PUD is similarly expected to be valid only for 
large or small values of /3i, the difference between the two regimes being only in the 
definition of the new timelike coupling 

A~^ , A large , (11) 



1 



(2A^2w*-i 



A ~ .o ".;a..,i ' A small . (12) 



In the large (3t limit the Migdal-Kadanoff approximation coincides with the ap- 
proximation scheme for high temperature QCD known in the literature as "complete 
dimensional reduction" . This consists in the assumption of complete decoupling of 
the non-static modes in the compactified time direction. The resulting theory is 
a (i-dimensional gauge theory coupled to a scalar field in the adjoint representa- 
tion which represents the fluctuations of the Polyakov loop around the minimum 
V{x) = 1. The corresponding action is the sum of the purely space-hke part of Sw 
and a new term Sh which is the remnant of the timelike part of Sw 

5, (mo) = § TrfE (-ml<f>{xy + Y.U,.,<f>{x)Uy{x + i)\ (13) 

X \ i=l / 

where 0(x) is an Hermitian N x N matrix and rriQ = d. 

This action is of the type studied by Kazakov and Migdal in the context of 



"induced QCD" [13| and in the present context by the authors of |14|. It is however 
well known that complete dimensional reduction does not take place in general and 
that non-static modes induce at one loop level some well defined interaction in the 
static sector. We expect this type of corrections to be the main source of error also 
in the Migdal-Kadanoff approximation, however it can be shown that they are of 
higher order in large N limit and therefore we shall not take them into account in 
the present paper. 

Let us conclude by noticing that the effect of the Migdal-Kadanoff approximation 
is to reduce the lattice size in the time direction Nt to its extreme value A'"t = 1. 
It will be important to keep in mind this fact in the last section, where we shall 
compare our prediction with the results of the Montecarlo simulations. 
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2.2 Decoupling of spacelike plaquettes 

An interesting feature of eq.® is that its timelike and spacelike sectors are only 
very weakly coupled. Let us try to make this statement more precise. First of 
all one can see that in the framework of a standard mean field approximation the 
decoupling of timelike and spacelike degrees of freedom in @ is complete f^. This 
result appears to be confirmed beyond the mean field approximation by montecarlo 



simulations 0. The decoupling was also checked in [^, where the contribution of 
timelike degrees of freedom to the space-like string tension at high temperature was 
evaluated in (2+1) dimensions within a strong coupling expansion and found to be 
neglected. This result turned out to be in agreement with high precision Montecarlo 



simulations 16 



In this paper we are addressing the problem of finding the critical deconfinement 
temperature and the relevant degrees of freedom appear only in the timelike part 
of the action. So, in agreement with the arguments given above, and following 
the literature on the subject [^, ^, |[, we shall from now on neglect the spacelike 
part of the action, eq.®. It is of course far from obvious to what extent such a 
crude approximation still holds at temperatures of the order of the deconfinement 
transition, so an independent check of its validity (such as results from Montecarlo 
simulations) is going to be important. 

The main outcome of this approximations is that in the resulting model the 
spacelike degrees of freedom can be explicitly integrated out, to obtain an effective 
action for the Polyakov loops only. This integration can be made (see section 5) 
by using a generalization of the Itzykson-Zuber-Harish-Chandra integral 0]. If 



one keeps only the lowest order in the coupling one recovers the action Seff, eq.(|TD 
discussed in the introduction. 

^e//(J) = E ^4^E Trf(\/(f))Trf(V^nx + z))| . (14) 

Here we have introduced the new coupling J defined in terms of /3 as: 

J = Pt/N' . (15) 

J has a smooth large N behaviour and so it is the natural coupling in the large 
A^ limit. In the following two sections we shall sistematically use J, referring to 
eg. ([TED for the connection with (3, and shall come back to the original coupling j3 
only in the last section, where the comparison with the Montecarlo results is made. 
Higher order corrections to eq. (14) are discussed in Section 5 where a similar 
analysis for the Heat Kernel action is also carried out. As we shall see, a phase 
diagram much richer than that of Seff emerges if one keeps the whole result of the 
integration. 



2.3 Large N Limit 

The model described by Sejf can be solved exactly in the large N limit 0], for any 
value of the space dimensions d, leading to a phase diagram with a first order phase 
transition located at J = 1/d. It is even possible to solve exactly the more general 
model in which the Polyakov loop is also coupled to an external "magnetic" field h, 

SeffUh) = E ^ej-^E Trf(K(f))TVf(Vt(f+z))| + 

+ /iA^^[Trf(\/(f) + V^(f))] . (16) 

X 

In this more general framework it becomes apparent that the first order phase 
transition is the end point of a line of third order phase transitions of the Gross- 
Witten type, located along the line 

I — 2h 
J = — -— , h,J>0 . (17) 

2.4 Comments 

There are two interesting observations which must be made at this point. First 
it is possible to show that the three steps of approximation described above com- 
mute, in the sense that the order in which they are made can be changed, but 
the same phase diagram is found at the end. For instance, if we neglect from the 
beginning the spacelike plaquettes, then the Migdal Kadanoff result eq. (|) can 
be obtained exactly within a strong coupling expansion (namely integrating all the 
spacelike links except the ones in the lowest slice). Similarly one could start from 
the beginning within the framework of a large A^ Eguchi-Kawai model, then make 
a Migdal Kadanoff approximation, and find again the same phase diagram (see |^ 
and this same observation in |^). This denotes an internal consistency of the set of 
approximations used which suggests that the value Jc = l/d for the deconfinement 
temperature should be taken seriously. Indeed, such value is rather stable under 
the various refinements of eq. (|l]) that will be discussed in the rest of this paper; at 
the end our results will never differ from it more than 20%. 

A second observation is that the magnetic field introduced in eq.(|l6D has a 
relevant physical meaning. It encodes the effect of the coupling of fermions to the 
pure gauge theory |18|, [1^. To be precise, one has to choose Wilson fermions and 



then make an expansion in the hopping parameter K. It is possible to see that at 
the lowest order in the K expansion the Zn symmetry of the pure gauge theory 
is broken and that the fermions contribute with a magnetic term h = 2N f{2K)^\ 
where Nf is the number of fermions introduced. 



3 The Heat Kernel Action 

The Heat Kernel action is best defined in terms of tlie cliaracter expansion 
exp(^H) = n 1 H drXr{Un-ij) exp — — |- \ X 

n,(i<j) V r \ '^Ps J } 

X n \J:drXriU^■,0^) eXp (-^) J (18) 



wliere tlie sum is over all the inequivalent, irreducible representations, labelled by 
r, dr is their dimension, Xr{U) is the character of U in the representation r and 
C^^^ is the value of the quadratic Casimir operator in the representation r. The 
index s and t label, as before, the spatial an temporal couplings respectively. The 
Heat Kernel action in eq.(|T8|) is normalized so as to coincide in the continuum limit 
{Ps, Pt -^ oo) with the Wilson action eq.(0). Hence in this limit eqs.(0J^J^) still hold 
with the substitution P — > P^ 

In the case of the Wilson action, both the Migdal Kadanoff approximation (sec- 
tion 2.1) and the decoupling and subsequent integration of the spatial degrees of 
freedom (sect. 2.2) require a truncation of the resulting actions so as to obtain the 
final form of the effective action eq.(|l^). The effect of these truncations can be 
neglected if Pt is very large or very small, but it is important in the intermediate 
region in which we expect the transition to occur. The situation is different for the 
Heat Kernel action. In fact due to the well known properties of the Heat Kernel 
action the strong coupling expansion which leads to eq.(|^) is exact in the timelike 
direction and the scaling of the Heat Kernel coupling is trivial: 

/5f -^^ = % (19) 

for all values of Pf^ and no truncation needed. As mentioned above, near the 
continuum limit (large P) the Heat Kernel and the Wilson action must coincide, 
in fact the same scaling law appears in eq. (|lT]) and (plQ]). On the contrary in the 
strong coupling regime the two behave very differently (as can be seen by comparing 
eq.(pJD and eq.(p!2|) ). This can be checked by a direct comparison of the strong 
coupling expansion of the Wilson and the Heat Kernel action (see section 5). The 
strong coupling expansion provides a relation between P^ and /5,which, as expected 
compensates exactly the difference between the two scaling laws eqs. (|19]), (|12D . 
This not trivial relation must be kept into account if one wants to make contact 
with say, results coming from Montecarlo simulations which are always made with 
the standard Wilson action. 

The decoupling of the spatial plaquettes, which we shall assume also in this 
context, makes it possible to integrate exactly on the spatial links by using standard 
properties of the characters. We have 

Zh = 



l[dV{x)YldU.{x) n It. drXr {u.{x)V{x + z)UhS)V\x)) exp (-^^) 



|l[dVix)l[J:XriVix + ^))xr{v\x))expL^^) (20) 

X x.i r \ / 



The effective action for tlie Polyakov loop given by (|20D , altliougli much more 
comphcated than that of eg. ([T6|) , can still be studied in the large N limit. This 
was done, both in the weak and in the strong coupling region, by K. Zarembo in a 
recent paper ||T0|. A rather puzzling result of [^ is the logarithmic dependence 
on the dimensions of the critical point, more precisely 

"^^ " 21og(2rf-l) ^^^^ 

This is a strong coupling result, obtained by looking at the point where the sym- 
metric vacuum becomes unstable, and it is apparently in disagreement with the 
results of the previous section which give a 1/d behaviour for the critical coupling 
Such disagreement is completely eliminated if one takes into account the relation 
between f]^ and (3. 

In the weak coupling region the action ( pOD can be studied by reducing it to 
a Kazakov-Migdal model with a quadratic potential, whose solution is given by 
a semicircular distribution of eigenvalues for the Polyakov loop. In ref. [^ it 
was remarked that an instability in the solution appears for the value of the cou- 
pling constant where the radius of the distribution develops an imaginary part, 
presumably a signal of some type of phase transition. Here we will produce a better 
understanding of this phenomenon, which will also allow a more precise determina- 
tion of the critical coupling. In order to do so we notice that the coupling between 
V{x) and V{x + i) at the r.h.s of eq. (^) is just the partition function of QCD2 on 
cylinder whose area Q ^ is ^73- and whose fixed holonomies at the edges are V{x) 
and V{x + i). 



In a recent paper [^ we have shown that QCD2 on a cylinder has a third order 
phase transition of the same type of that found by M.Douglas and V.Kazakov in 
the spherical geometry. Such a transition can be seen in the language of 2d QCD as 
an istanton driven delocalization of the trivial semicircle distribution of eigenvalues. 
We will show in the next section that in the present context of Finite Temperature 
LGT this phase transition is responsible for the instability of the weak coupling 
solution. Before proceeding to do that however, we shall end the present section 
by reviewing and completing the results of ref. [^ on the Douglas-Kazakov phase 
transition on a cylinder. 



^As usual in the contest of QCD2, the area A of the cyhnder is expressed in units of \, where 
§2 is the QCD2 coupHng. 
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3.1 Douglas-Kazakov phase transition on a cylinder 

As previously discussed the Heat Kernel effective theory for the Polyakov loop is 
given in terms of the partition function /C2 (fi'i , ^72 ; -^) = ^^2(0,^;-^) of QCD2 on a 
cylinder, with fixed holonomies (Polyakov loops) gi and g2 (whose invariant angles 
are {0j} and {9i} respectively) at the two boundaries. 

This partition function is the solution of the Heat Kernel equation on the SU{N) 
group manifold : 

[n^ - \j-\<\>) E ^J{<\>) - ^Ar(iV^ - 1)) /C,(0, e,A) = ^ (22) 

uniquely determined by the condition : 

hm ;C2(5'i,fi'2;^) = K9\.92^) (23) 



where, in the notations of ref. |^, i7(0) is the Vandermonde determinant for a 
unitary matrix: 

J(0) = n2sin^i^ (24) 

and b is the invariant delta function on the group manifold. An explicit form of 
/C2(0, Q\ ■^ is usually given by the standard character expansion : 



)C2{g\,g2\A) = E^xp 






2N 



xMXriO) (25) 



As pointed out in [^, ^, upon redefinition of the Kernel by /C2(0, 9', A) ^ }C2 = 
i7(0)i7(6')/C2, eq. (p^) becomes the (euclidean) free Schrodinger equation for A^ 
fermions on the circle, where A plays the role of the (euclidean) time. This means 
that, because of the condition (|23|), we can interpret lC2{(p,0]A) as the euchdean 
transition amplitude of this system of fermions from the configuration {0j} at the 
zero time to the configuration {—9i\ at the time A. 



By using the results of ref. [53] and [21| the above expression can be rewritten 



as a sum over the A^ integers U labelling the co-root lattice of SU{N): 

/7V\i/2 /A \ (A){i-N)/2 



/r.\ , iV(iV-l) 



l)-(m^T-^^exp 

ft} 



2^.=i 



N 



(26) 



where P denotes a permutation (of sign cr(P)) of the indices. The meaning of the 
modular transformation relating this last expression to the character expansion is 
clear: the integers labelling the unitary representations in the character expansion 

11 



(p5|) correspond to discretized momenta of the fermions on the circle, while eq. ( p6| ) 
gives the corresponding coordinate representation and the integers in the co-root 
lattice are interpreted as the winding numbers of the fermions. 

The large A^ behaviour of the partition function (|25|) in terms of the eigenvalues 
distributions po{x) and pi{x) corresponding respectively to {(pi} and {9i} can be 
obtained starting from the differential equation (^) and then applying the same 
procedure used by Matytsin in |^ to find the large N limit of the Itzykson-Zuber 
integral. The time (area) evolution of the eigenvalue distribution p{x) is given by 
a Das-Jevicki hamiltonian 



25 



H[p{x),n{x) 



dxp{x 




TT" 



p'ix) 



(27) 



where n(x) is the canonical momentum conjugate to p{x). The only difference with 
respect to ref. [^ is that the density p{x) is here a periodic function of period 27r. 
In terms of the complex quantity 



f{x,t) 



dU{x,t) 
dx 



+ inp(x,t) 



(28) 



the equations of motion become the Hopf equation of motion for an ideal fluid : 



dt dx 



0, 



with the boundary conditions 

p{x,t = 0) = npo{x) , p{x,t = A) = npi{x) 
A solution to eq. (l29| ) can be obtained from the ansatz 



(29) 



(30) 



p{x,t) = 
p{x,t)=0 



nr^it) 



r\t) 



x^ 



\x\ < r^{t) 
r'^{t) < \x\ < n 



(31) 



The ansatz ( pT]) provides a solution of the Hopf equation if the time dependence 
of the radius r of the distribution has the form 



r(t) 



'(t + «)(/?-t) 



a + P 



(32) 



The arbitrary constants a, (3 are determined by the boundary conditions (pOl). Nat- 
urally the initial and final distributions in (|30| ) have consistently to be semicircular 
with radia r(0) = ro,r{A) = ri. Given tq and ri, the radius of the distribution 
is determined at any section of the cylinder by eq. (P^) . However because of the 
periodicity condition on the eigenvalue distribution the density ( |3T1 , |32D is a solution 
of the saddle-point equations (2£,^) only if r(t) < n for any t on the trajectory. 
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For any given value of tq and ri, the maximum value of r{t) increases as the area 
A increases , so the solution ( ^|32D is valid only if the area A is smaller than a 
critical value Ac, where the maximum radius equals vr and the eigenvalues fill up 
the whole circle. The critical value of A at which such phase transition occurs can 
be easily calculated and is given by 



^c 



\ 



7^4 _ 7^2 



rg + r^ 



- + y (tt^ — Tr^rori)^ — tc^Itq — ri)^ 



(33) 



Note that the partition function of QCD on a sphere is a particular case correspond- 



ing to tq = ri = 0. Indeed in this case the critical value 



becomes Ar 



TC 



which is just the value found by Douglas and Kazakov ||TT|. Hence the large A^ 
phase transition between the small area (gaussian) phase and the large area phase 
is the generalization to the cyilinder of the Douglas-Kazakov phase transition on 
the sphere . 

The meaning of the phase transition was discussed in pO[ : in the gaussian phase 
the eigenvalues are confined in the (— vr, vr) interval, namely the configurations where 
a fraction of the eigenvalues wind around their configuration space do not contribute 
in the large A^ limit. Correspondingly all integers k in eq. (^) can be set to zero. 
Such topologically non trivial configurations contribute instead in the large N limit 
of the large area phase. 

The role of instantons in inducing the Douglas-Kazakov phase transition on the 
sphere was further investigated by Gross and Matytsin p6|. 



4 Weak Coupling 

As a result of the approximation schemes discussed in the previous sections, we 
can now write an effective action for the Polyakov loop starting either from the 
Wilson action or from the Heat Kernel action. Moreover in the large A^ limit, 
we can assume that the saddle point solution is translational invariant and, as a 
consequence, the action reduces to a one plaquette integral, both in the Wilson and 
in the Heat Kernel case. The partition functions in the two cases are given by 



i=l 



dV 



NJ Re Tr UVVW^ 



and 



'H 



dV 



dUe 



}C2{V,V^;1/J'') 



(34) 



(35) 



where d is the number of space dimension and /C2(V^5 V^; 1/J^) is the Heat Kernel 
on a cylinder discussed in the previous section. The couplings J and J^ are related 
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to the coupling f3t of the theory defined on a lattice with Nt links in the time 
direction by eqs. (|TT| , |12|, |l^) and keeping into account the redefinition (^). We 
shall study the phase diagrams resulting from (^^ and (|35|) keeping in mind that, 
in the region where the two actions describe the same physical system, the same 
critical phenomena should occur at values of J and J^ corresponding to the same 
value of Pt- On the other hand, we know that the rescaling of J^ given by eq.(|19D 
is exact after decoupling the spatial plaquettes. Thus the comparison between the 
two partition functions can be used to improve eqs (0) and ([T^). In the case of 
eq.(|TTD, for instance, we will find that the two J's must differ by a constant term 
for the two theories to match in the weak coupling region. 



4.1 Weak coupling expansion for the Heat Kernel action 



The partition function (|35D was studied in detail by Zarembo in |jTO[, and part of 
our analysis will overlap with his work. We will, however, be able to interpret the 
critical phenomena occurring in the weak coupling region in terms of the Douglas- 
Kazakov phase transition in the underlying QCD2 structure. Let us first rewrite 
eq. {EE) explicitly, as an integral over the invariant angles of V. 



'H 



X{deAj\e) 



(i-d) 



k, P 



.IY(P)(^-\NJH ^^(e,-ep„+27rfc,)2 



(36) 



As discussed in the previous section, if J^ is larger than a critical value (to be 
determined), the contributions of the winding configurations can be neglected in 
the large N limit, and all fcj's in eq.(|5BD can be set to zero. Moreover, for small 9, 
we have 

j2(^) = A2(0)e-ii^S«^?+«(''') . 

Inserting these results in eq.(|36|) we obtain 

r id 



(37) 



n^^* A'w 



(i-d) 



y-i^^iyin^^Nj^e^ept^i) 



(3J 



Provided J^ is larger than its critical value, the only approximation needed to 
go from eq.(^) to eq.(^) consists in neglecting 0{9'^) terms in (^7|). The r.h.s. 
of eq.(^) is a Kazakov-Migdal model with quadratic potential in the mean field 
approximation, which was solved exactly in the large N limit, for any number of 
space dimension, by Gross P^. The eigenvalue distribution is a standard Wigner 
distribution 

(39) 



p{d) 



yr2 - e^ 



nr^ 



where the radius r of the distribution is given by 
2 4(2rf-l) 



J"{m?{d - 1) + dJm^ - A{2d - 1)) 



(40) 
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The parameter m which appears in ( ^0] ) is the usual mass term of the Kazakov- 
Migdal modeL Following the standard normalization we have that the coefficient 
of the quadratic term in 9 in eq.(|38D must be ^pr. This implies 



m 



2d 



d 



QJ^ 



(41) 



In 



rO[ , Zarembo argues that a phase transition occurs at a value of J^ equal to 
(or larger than) the one for which r^ becomes complex. Up to the effect of terms 
of higher order in ^, we can in fact determine the exact value of J^ at which the 
phase transition occurs, and we will further argue that the phase transition is of 
third order. From the analysis of the previous section, we know that the kernel 
^2iV, V"^; l/J^) undergoes a phase transition of the Douglas-Kazakov type when 



J^ 



4 2 2 



(42) 



where r^ is given by (0). The results for the critical value of J^ , and for the value 
Tc of the radius at the critical point, are given for different space dimensions in Table 
(I). The critical values of J^ obtained by Zarembo are also listed for comparison. 



d 


re 


J^, DK 


J^'.z 


2 


2.96 


0.321 


0.311 


3 


2.80 


0.226 


0.218 


4 


2.66 


0.192 


0.184 


CX3 


0. 


I/tt^ ~ 0.101 


1/12-0.083 



Tab. I. Value of J^ below which the Wigner distribution becomes unstable, from 
the Douglas-Kazakov phase transition (third column) and from ref / |7^/ . The 
radius of the distribution at the critical point is given in the second column. 



We observe that our values are consistently just slightly higher than Zarembo's, 
which clearly indicates that we are dealing with the same instability and not with 
two different physical phenomena, and on the other hand shows that the present one 
is a better approximation. As a consistency check, we notice that the critical radius 
is less than vr for any number of dimensions, whereas if we take the solution of [^ 
we see that it exceeds tt for d = 2 and d = 3. In spite of the small value of J^ we 
expect the weak coupling solution to be very reliable above the Douglas-Kazakov 
phase transition particularly for large d, where the small value of Vc ensures that 
the quartic terms which have been neglected in eq. ( p8| ) are indeed small. 

To sum up, we have established that the Wigner distribution of eigenvalues be- 
comes unstable at a value of the coupling constant given by (|42|), and the transition 
is driven by the winding modes in the configuration space of the eigenvalues. This 
analogy with the Douglas-Kazakov phase transition on a sphere gives an almost 
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compelling argument that the corresponding phase transition is of the third order. 
The appearance at the critical point of new classical trajectories corresponding to 
winding eigenvalues has presumably the effect of spreading the distribution, espe- 
cially at the extremes. On the other hand, for low dimensions, this transition occurs 
at a radius of the eigenvalue distribution very near to tt. This probably means that 
the maximum corresponding to the broken, deconfined phase becomes unstable, and 
the distribution would collapse into the uniform one. However this does not mean 
that the Douglas-Kazakov transition can be identified with the deconfinement tran- 
sition. This is not the case for high dimensions, where the critical radius is small, 
approaching zero as d increases. In this case the Douglas-Kazakov phase transition 
is a transition from a classical Wigner distribution to another one, so far unknown, 
but still presumably peaked around the origin. 



4.2 Weak coupling expansion for the Wilson action 



The partition function Ziy of eq. (J3^) can be written as an integral over the invariant 
angles of the Polyakov loop V, 



Zw = JTldea\e){x{e,j)r 



(43) 



where 



X(0,J) 



r TTji AfJ Re Tt {U cos 9U'' cos e+U sineUUm9) 



(44) 



and we denote by cos 9 and sin 9 diagonal matrices with eigenvalues cos 9i and sin 9i 
respectively. In the large A^ limit this integral can be evaluated as a sum of saddle 



point contributions ||28[. The result is 



^0, J) = E 



exp 



NJ J2i cos( 



^P{i)) 



Hp{9) 



^^+^W) 



(45) 



where 



Hp{9) = Y\_ (cos 9i — cos ^j) (cos 9p{i) — cos 9p(^j)) 



i<] 



+ (sin 9i — sin 9j) (sin ^p(j) — sin 9p(^j)] 



(46) 



The denominator Hp{9) in eq. (|45|) has a double pole whenever two ^'s coincide, 
so it can be written as A^(6') exp {hp{9)), and hp{9) can be easily calculated for small 
9, with the result 



hpi9) = -f E ^' + ^ E ^^^P» + Oi9'') 



(47) 
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For small 6, the partition function can thus be written as 

id 



'W 



-' i I p 



(48) 
Once again, this is a Kazakov-Migdal model with quadratic potential, and it coin- 
cides with (eq. (|38D ) if we identify 

J = j'' + \ . (49) 

We have thus identified the leading correction to eq. (0) away from the con- 
tinuum (weak coupling) limit: the two lattice theories defined by the Wilson and 
the Heat Kernel actions begin to differ along the respective renormalization group 
trajectories by a constant shift of their couplings. As J decreases higher order terms 
in the angles 9 will become relevant and the two theories will become increasingly 
different. 

As for the phase diagram, one cannot apply to the Wilson action the same pre- 
cise estimate based on the Douglas-Kazakov phase transition, but only the simpler 
argument given by Zarembo. This is due to the fact that the Douglas-Kazakov 
transition crucially depends on the precise form of the Heat Kernel action, and in 
particular on the winding numbers ki which are neglected in the Kazakov-Migdal 
approximation. In any case, as mentioned above, the phase diagram of the two 
theories is the same as long as the relation (ESI) is duly taken into account. 



5 Strong coupling 

While the inverse coupling jSt/N"^ in the lattice theory with A''^ links in the time 
direction becomes very large as we approach the continuum limit, the effective 
coupling J in the reduced A^ = 1 theory ( ^4]) remains relatively small, typically 
of order 1/d, even at temperatures approaching the deconfining transition. An 
expansion in powers of J is therefore appropriate to study the critical phenomena 
in the effective models for the Polyakov loop considered in the previous section. 

As a preliminary step let us establish some notations. We are interested in the 
large A^ limit, so the fundamental quantity is the distribution p{6) of the eigenvalues 
of the Polyakov loop. It is convenient to expand p{6) in its Fourier modes 

-| OO 1 oo 

Pid) = ^ E Pne'"-' = ^ E xne'---'^-' , (50) 

n=— oo n=— oo 

where a„ G (— 7r/2,7r/2) is the argument of pn modulo vr, and x„ coincides with 
the modulus of pn up to a sign 0. The reality of p{9) requires p_„ = p* , while the 
normalization of p{6) to 1 in the interval (— vr, vr) fixes po = 1- 

^It is convenient for the following discussion to restrict the range for a„ and have Xn taking 
also negative values 
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The inverse formula 



Pr, 



p{e)e 



-ind 



(51) 



shows that p„ corresponds to the large N limit of the loop winding n times in 
the time-like direction with the given eigenvalue distribution. In particular, p±i 
corresponds to the large N limit of the Polyakov loop. 

The Zjv invariance of the effective theory becomes, in the large A^ limit, a U{1) 
invariance under the shift 9 ^>- 9 + 5, that is «„—;►«„ + n5. If this symmetry is 
unbroken the eigenvalue distribution is simply given by p{9) = -^^ namely x„ = 
for n 7^ 0. In the broken phase the f/(l) symmetry connects the different vacua. 
If we choose the vacuum peaked at = 0, then the symmetry of the action for 
9i -^ —9i will force the vacuum distribution p{9) to be even in 9, thus fixing all «„ 
to zero. In this situation the eigenvalue distribution takes the form 



"f^' - T. 



l + ^Xn 



COS nc 



n=l 



(52) 



with Xn real. 

The next ingredient is the integration measure J'^(9\ which in the large A^ limit 
can be expressed as a function of the loop variables p„, as 



J'' {9) = exp 



lim -7^ d9 dip p{9) p{ip) log(l - y cos{9 - if)) 



(53) 



where the double integral, which would be divergent a.t 9 = ip, has been regularized 
by the inclusion of the parameter y, and terms in the exponent suppressed by powers 
of A^ have been neglected. Eq. ( ^3]) can be calculated by expanding in powers of 9 
and p and resumming the resulting expression. The result is 



J^i9) = exp 



\unN'ICo{y) + Y.Ck{y)xl 



k=l 



(54) 



where Co{y) is an irrelevant divergent expression and Ck{y) is given by 



Ck{y) 



y 



(55) 



After removing the divergence, the limit y ^ 1 can be taken, and gives 



J' (9) = exp 



N" E u^k 



fc=i 



k 



(56) 



Let us consider now the partition function Zw, given in eq. (^31) . We are going 
to show that the integral T{9, J) can be written in the large A^ limit as 



j(0,j) = e^'i:r.i^^-^^(-.-) 



(57) 
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where the functions Fk{x, a) can be shown to depend only on Xi and at with i < k, 
and can be determined in principle by using Schwinger-Dyson equations. 
Correspondingly, the large N limit of the partition function Zw becomes 

1 °° / t2\ 

— In Zw = j:[d J' F,ix, a)- ^j . (58) 



It is already obvious from ( pq ) that, for small enough J, the free energy has a 
maximum for Xk = 0, that is for a uniform distribution of eigenvalues. The first 
order approximation of ( |58D is well known (see for example [Q), and is simply 
Fi{x, a) = x\. The free energy to this order becomes 

-^\nZw^{Jd-l)x\ , (59) 

with all other Xk set to zero. The phase structure of (|59D is very simple: for J < ^ 
the maximum of the free energy occurs at Xi = and we are in the unbroken phase. 
Above the critical point J = ^ we have instead Xi = 1/2, which is the maximum 
value allowed by the positivity condition on p{6). This is exactly the phase diagram 
of the large N limit of the effective action (|I]) discussed at the beginning of the 
introduction. In fact it is easy to see that at this first order approximation the two 
actions are actually equivalent. This allows to better understand in which sense we 
expect, going to higher order in eq. (|58D , to improve the effective action (0). 

Let us compare the first order expansion of the Wilson action with the corre- 
sponding one for the Heat Kernel action Zhk, given in eq. ( p5|) and (|36|) . The 

strong coupling expansion for the kernel on the cylinder coincides with the charac- 

1 
ter expansion, which is a power series in e 2jh . The truncation of the series to the 

first order corresponds to a truncation of the character expansion to fundamental 

representation, namely 

c(2) 1 

J2xrio)xri-e)e-^^ ^ i + 2xfie)xfi-e)e-^ 



exp 



N^2e~^J^xl , (60) 



where in the last step we have used the fact that X/(^) = Nxie^°'^. At first order 
in the strong coupling expansion Z^ and Z^k coincide, provided we make the 
identification 

J = 2e'^ , (61) 

in agreement with eq.s (|11|) and (|12|). Indeed, we have just rederived, comparing 
the first order strong coupling expansions of Zw and Zhk, the well known rescaling 
properties of the coupling constant J. 
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5.1 Strong coupling expansion for the Wilson action 

The functions Fk{x,a), which characterize the strong couphng expansion of Zw, 
can be determined by solving Schwinger- Dyson type of equations. Let us consider 
the integral 

J( J; $, vj/) = |[rff/]e^^^K^*i^t*i+^*2^t*2) , (62) 

where $„ and \&q are hermitian matrices. This integral is more general than the 
one in eq. (0), which is obtained from X(J; $, \E') by imposing the constraints 

E*^ = E^^ = 1 (63) 

and by identifying $ with ^. Notice that, for $2 = ^2 = 0, eq. (|62D reduces to the 
well known Itzykson-Zuber integral. The Schwinger- Dyson equations for X( J; $, \E') 
can easily be proved; they read 



Tr 



did Id 



Nd^^.Nd^^, Nd^^^ 



X( J; $, ^) = rrr (vl/„^vl>„^...vl/„J X( J; $, v]/) 



(64) 
These equation could, in principle, be solved order by order in J. It is however 
more convenient to make use of the fact that X( J; $, \1') reduces for $2 = ^2 = to 
the Itzykson-Zuber integral, whose "strong coupling" expansion has been explicitly 
derived in |jT^ up to the eighth order in J. In fact it is not difficult to verify that. 



given a term in the expansion of the Itzykson-Zuber integral, containing only $1 and 
^1, there is a unique extension of it to include the dependence on $2 and \E'2 which 
satisfies the following requirements: i) it is invariant under the 0(2) symmetry of the 
integral ( |62D and of the constraints (p3|); ii) it contains in each term $j and \l/i with 
the same power, and thus is invariant when $j and \E'i are rescaled respectively by a 
factor A and j, which is a symmetry of the integral before imposing the constraints 
(|6^). Using this property, it is easy to generalize the known results concerning the 
Itzykson-Zuber integral, and calculate the functions Fk{x, a) appearing in eq. (p7\). 
Up to fc = 4 they are given by 



Fi{x,a) = Xi , (65) 

F2{x,a) = - n — 2a;^ + 2xf + X2 — 2x^X2 cos(2q;i — 0:2)) , (66) 

F'i{x,a) = —fSx^ — 12x^ + 16x^ + 12x^^2 + X3 + 6x^x2 cos(2ai — ^2) 
—24x^X2 cos(2ai — ^2) + 4x'^X3 cos(3ai — a^) 
—6x1X2X3 cos(a;i + 0:2 — as) ) , (67) 
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F4{x, a) = ^ (-1 + 34x^ - 144a;? + 192x1 - mx\x\ + 226xtx^ 

+4^2 + 20x1X3 + X4 + 34x1X2 cos(4ai - 202) 

— 12x^X2 cos(2ai — ^2) + 152x^X2 cos(2ai — 02) 
— 384xiX2 cos(2ai — 0^2) — 48xiX2 cos(2a;i — 0:2) 

— 16X1X3 cos(3ai — 0:3) + 8OX1X3 cos(3a;i — 0:3) 
+8x1X2X3 cos(ai + ^2 — as) — 120x^x2X3 cos(q;i + 02 — as) 
+16x1X2X3 cos(ai — 2^2 + as) — 10X1X4 cos(4ai — 0:4) 
+20x1X2X4 cos(2ai + 0:2 — 04) — 4x2X4 cos(2a2 — ^4) 

— 8x1X3X4 cos(«i + ^3 — a4)) . (6J 



Notice the invariance of the functions -Ffc(x, a) under a^ ^ a^ + ^5, as a result 
of the f/(l) invariance of the theory in the large A^ limit. 

We can now substitute the expressions (|65| - ^8[) into (H), and study the free 
energy up to fourth order in J. We have already noticed that at first order a phase 
transition occurs at J = i. The higher order terms in J can be regarded then as 
higher order terms in i, and the strong coupling expansion as a large d expansion. 
Consider for example the coefficient of the quadratic term in xi. The vanishing of 
this coefficient signals an instability of the symmetric vacuum x^ = 0. This leads 
us to study the equation 

ciJ-^j2 + ^J='-l + O(j5) = 0, (69) 

which can be solved order by order in i, leading to 



It is instructive to compare this result with the result obtained by Zarembo 
using the Heat Kernel action. He finds that the exact value at which the instability 
occurs is given by J^ = 210 (L-i) ' When inserted in ([6T| ) this expression gives 
J = ^_\/2 ^ whose expansion in powers of 1/d coincides with the one found in ([70| ) 
up to the third order. This is more that what one would have naively expected, as 



eq. (Ul) was derived by equating the expansions of the Wilson and the Heat Kernel 
actions only at first order. 

While eq. ( [70| ) gives an upper limit for the critical value of J, finding the 



precise value at which the deconfinement transition occurs, requires a more detailed 
analysis. As discussed at the beginning of this section, when searching for the broken 
vacuum we can set all a„ to zero, and regard the free energy as a function of the 



Xn only. By using eq. (^ - |6g) we obtain 



1 dJ'^ 

In Zw = dJxl + -— (1 - 2xi + 2xi + X2 - 2xiX2 



Ar2 "^ -^ ' 4 
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+ — - (Sxl - 12xt + IQxl + V2x\x\ + x\ 

+ 6x^X2 — 24x^X2 + 4x'^X3 — 6x1X2X3 

+ — - (-1 + 34x^ - 144x^ + 192x^ - 40x?x^ 

+ 226x^x^ + 4x^ + 20x?x^ + x^ + 34x^x^ 

- 12x1X2 + 152xiX2 - 384x^X2 - 48xiX^ - 16x^X3 
+ 80x^X3 + 8x1X2X3 — 120x1X2X3 + 16x1X3X3 

— 10x1X4 + 20x1X2X4 — 4x2X4 — 8x1X3X4 

_ 2__2__2__2 /yi \ 

■^1 ^^2 0^3 7-^4 • V^) 



We must now look for the inaxiinum of (0), within the domain where p(^) is 
positive or zero for any Q. The equations for such a domain can be obtained in 
principle by requiring that, for some 6'o, 

p(^o) = ^pW|.=.o=0 (72) 

and by eliminating Qq from the equations. The resulting equations for the x„'s give 
the boundaries of the physical region. In order to have an intuitive picture of how 
the deconfining transition takes place let us set X3 = X4 = in (^) and plot the 
free energy as a function of xi and X2 at various values of the coupling J. This 
is shown in Fig.l, where a contour plot of the free energy is shown for d = 2 and 
J = 0.55, 0.60, 0.66. 

The physical region in the (xi,X2) plane can be easily determined, and it is 
represented by the region inside the thick line. The straight edge on the left is 
given by the equation xi — X2 = 1/2 and corresponds to density distributions 
vanishing at 6* = vr, whereas the the straight edge on the opposite side corresponds to 
distributions vanishing at ^ = 0. The plot at J = 0.55 clearly shows the maximum 
at Xi = X2 = 0: the system is in the unbroken phase dominated by a constant 
distribution of eigenvalues. In the next plot, at J = 0.6, a local maximum has 
appeared at the edge of the physical region, and it is becoming competitive with 
the unbroken maximum. Notice that there is a symmetry xi -^ — Xi, so that 
a symmetric maximum appears on the other edge of the physical region. This 
symmetry is accidental, and it is removed when X3 and X4 are switched on. In the 
last plot, at J = 0.66, the maximum at Xi = X2 = has disappeared and the system 
is clearly in the broken phase. 

When the dependence of (^Tj) on X3 and X4 is taken into account, the boundary 
of the physical region is given by fourth order algebraic equations, together with 
the hyperplanes xi — X2 + X3 — X4 = 1/2 and xi + X2 + X3 + X4 = —1/2. The phase 
transition occurs when the maximum value of the free energy on the boundary 
becomes larger than the value at the symmetric point. The critical point can be 
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Fig.l Contour plots of the free energy for: (a), J = 0.55; (b),J = 0.61; (c), 
J = 0.66. The numbers 1 and 2 mark respectively the presence of the symmet- 
ric maximum (xi = X2 = 0^ and the "broken" one. The thick line encloses the 
"physical" region 
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determined numerically, and it is given for different dimensions by the values listed 
in Table II. 



d 


Jc 


Xi 


X2 


X3 


X4 


2 


0.601 


0.50 


0.06 


-0.03 


-0.03 


3 


0.379 


0.52 


0.03 


-0.03 


-0.02 


4 


0.275 


0.53 


0.02 


-0.03 


-0.02 


100 


0.010 


0.49 


0.004 


0.004 


0.005 



Tab. II. Values of the critical coupling Jc at the fourth order in the strong 

coupling expansion of the Wilson action. In the last four columns the 

corresponding values of Xi, [i = 1 — A) are reported. 

The values of Jc depends of course on the order at which the strong coupling 
expansion is truncated, the lowest order being, as already mentioned , 1/d. The 
second and third order give for d = 2 J = 0.66 and J = 0.62 respectively, showing 
that the value reported in Tab. II is likely to be an upper bound. Higher dimensions 
are much less sensitive, in agreement with the idea that strong coupling expansions 
is also a large d expansion. 

It would be interesting then to discuss these results also from the point of view of 
the 1/d expansion, along the lines followed to study the instability of the symmetric 
vacuum. However, as the vacuum corresponding to the broken phase occurs on the 
boundary of the physical region, the constraints (|72D must be taken into account, 
giving an extra effective potential of order 1 in 1/d. This means that already in 
the zeroth order of the expansion the Fourier components Xn with n > 1 are not 
zero in general, and Xi does not coincide with its naive asymptotic value xi = 0.5 , 
as clearly indicated by Table II. The zero order approximation is already far from 
trivial when dealing with the broken maximum near the boundary, and we will not 
attempt its solution here. 

5.2 Strong coupling expansion for the Heat Kernel action 

We consider now the problem of obtaining a strong-coupling expansion for the 
effective theory of the Polyakov loops using the Heat Kernel action instead of the 
Wilson action. 

The starting point is naturally the character expansion 



Cr 



JC2{gu92'; ^/J") = T.^'^^XR{9i)xR{g2^] 

R 



(73) 



Considering the translational invariance of the saddle point solution for the effective 
theory, we are actually concerned with the case gi = g2 = V, Trl^ being the 
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Polyakov loop. However, we shall consider the more general case gi ^ (72 and put 
91 = 92 = V only at the very end. 

The large A^ limit of the kernel (^) was studied in Section 4 for large J^ . The 
"gaussian" solution which arises in that situation is valid only down to a critical 
value of J^ , where a third order Douglas-Kazakov phase transition takes place. It 
would be interesting to find out the exact expression in the large N limit of the 
kernel ([73|) , also in the strong coupling phase, extending the results of ref. [|lll] for 
the sphere to the cylinder. 

Unfortunately, such an exact expression is not yet available. We shall instead 
give an expansion of the free energy F{gi, (7^^; J^) corresponding to the kernel ([73|) 
in powers of the relevant coupling, which in this case must be taken to be 



rJ qn ^O 



(74) 



(notice that at the first order J^ coincides with J). Define the free energy in the 
large N limit by 

K.i9i,92';jw) = e''''^'^^'^"'"^ , (75) 

where as usual terms suppressed by powers of A^ in the exponent will be neglected. 
We can write then 



F{g„ g,\ J^) = 2 e-^'^-^" F^ {g,,g-\ J^) + 4 e-'/'"F,^{g„ g,\ J^) + . 



(76) 



where the residual dependence on J^ in the individual terms is, as we shall see, 
polynomial in 1/J^. 

The expansion (^) is of the same type of the one obtained when using the 
Wilson action; it is therefore what we need in order to compare the two cases. At 
lowest order, in fact, the two expansions coincide after the identification ([7^). 

We leave to the Appendix some technical remarks about the way to obtain such 
a strong coupling expansion, and about its meaning and its potential connection 
with a string representation for the gauge theory. Presently, we quote the result of 
the expansion of the free energy F(y, V^; J^) up to fourth order in Jyj. The result 
is expressed in terms of the Fourier modes of the density of eigenvalues p{9) , that 
is of the quantities Xn and a^, introduced in the previous subsection. We find, in 
the notations of the previous section. 



Ff (x, a; J^) 
F^ix.a-J'') 



F^ix.a-J'') 



1 

4 
1 
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—2x\ H ^^1 + ^2 ff7-^i'^2 cos(2a;i — 012) 



2 



(77) 
(78) 



Qxi 
12 



J^ 



2'^l 



JH 



4^1 + 



12 



Uj-lihr} "l ^iho 



+ — ttX'jxS cos(3ai — a^) 7732^12^2 cos(2q;i — ^2) 



J" 
6 



J^ 



rX'^xcos(3ai — 03) 



jH3 

12 
-jj^XiX2Xs cos(ai + 02 - «3) 



(79) 
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Ff(x,a;J^) 



1 
16 



-2x1 



10 



2 2 I 
2-^l"''2 + 



10 



X 1 



2J"^ 



x\ + 



JH 



4 4 I 
^XiX2 + 



4 
■—j^x\x2 cos(2ai — 02) 

4 
ra;2X2 cos(2ai — 0:2 
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4 + 



37^*^ 
6 



X 1 



Xn 



J^ 



2 2 I 

2"''l"''3 + 



6 



3J^' 

8 

9 



rX^X2 COs(2q;i 



Z'^A 



«2j 



X]^X2 cos(2a;i — 0:2) 



X1X3 cos(3ai — 0:3) H 774a;ia;3 cos(3ai — 03) 



2J^^ 



J^ 



X1X2X3 cos(a;i + 0:2 — 0:3) 



18 3 , 

—3X^X2X3 cos(q;i + ^2 - «3) 



6 



JH- 



iXiX^Xs cos(ai — 2^2 + c^s) 

;x\x2Xi COs(2ai + ^2 — «4) 



(■x^X4 cos(4ai 



— ^^1X3X4 cos(a;i + 0:3 — a4j 



3J^' 

—7^X2X4 cos(2a;2 — o.^ 



Oii) 



m) 



Some remarks are in order. First, as we mentioned above, each coefficient in 
the expansion of the free energy in powers of J^ still contains powers of jjt, that is 
of log Jw This is to be contrasted with the behaviour of the free energy associated 
with the Wilson action, where each term in the expansion is a function of x„ and a„ 
only. Second, the expansion ( [75| ) agrees to all available orders with the result of |TD|, 
that the symmetric vacuum becomes unstable at J 



H 



or J,, 



1 



21og(2d-l)' ^^ ""w d-lj2- 

2 \„ ''^), and following the 



This can be checked by looking at the coefficients of x{ in (^ 
arguments leading to ([70|) in the previous subsection. 

While the lowest order in J^ of the expansion ( pO]) coincides with the corre- 
sponding one for the Wilson action, the higher orders are different, although the 
general structure remains the same. It is interesting to see how the differences of 
the two models reflect in the critical values of J^, at various dimensions. The same 
type of numerical analysis that lead to Table II using the Wilson action now gives 
the results listed in Table III. The same critical values, in terms of the original Heat 
Kernel coupling J^ are listed in the fourth column of Tab. IV. 



d 


■Jw 


Xi 


X2 


X3 


X4 


2 


0.601 


0.41 


-0.13 


0.03 


0.07 


3 


0.339 


0.45 


-0.10 


-0.01 


0.05 


4 


0.238 


0.50 


-0.07 


-.05 


0.02 


100 


0.0096 


0.64 


0.07 


-0.10 


-0.03 



Tab. III. Same as Tab. II, but for the Heat Kernel action. Notice the coupling J^ 



is related to the Heat Kernel coupling through eq. / fZ^j 
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The coincidence of the critical value of Jw for d = 2 with the one given in Table 
II is probably accidental. However, it can be seen from eq. (|80D that some relevant 
coefficients in the first three terms of the expansion coincide with the corresponding 
ones for the Wilson action when J^ = 0.5, that is when J^, is around 0.75. This 
leads to expect a better agreement between the two theories in that region of J^. 
In dimensions larger than 2 the discrepancies are of the order of 10 — 15%, and they 
probably give an idea of the dependence of the results on the chosen regular izat ion. 



5.3 Summary of the results 

The results of both the weak and the strong coupling expansions are summarized 
in Tab.IV. 



d 




Jdk 




Jc 


^ S.V.l. 


J s.v.i. 


2 


0.321 


0.561 


0.416 


0.601 


0.455 


0.667 


3 


0.226 


0.468 


0.282 


0.379 


0.310 


0.400 


4 


0.192 


0.434 


0.235 


0.275 


0.257 


0.286 


oo 


4 -0.101 


0.333 


1 


1 
d 


1 


1 
d 


21og(2a!) 


21og(2a!) 



Tab.IV. Summary of our results. D.K refers to the Douglas- Kazakov point, 

s.v.i. to the symmetric vacuum instability point and Jc, J^ are our best estimates 

for the deconfinement phase transition. See text for detailed explanation. 

In the second and third columns the weak coupling results are listed, namely 
the value of J^ and J for which the semicircular distribution becomes unstable (as 
discussed in sect. 4), in the Heat Kernel and Wilson model respectively. 

In the two middle columns the best estimate for the phase transition, using 
the fourth order strong coupling expansion, is given for both models. In the last 
two columns the values of Jh and J at which the symmetric vacuum becomes 
unstable are given. They are obtained from the equations J^ = l/(21og(2(i — 1)) 
and J = l/{d — 1/2), the former is an exact result |]IU|, the latter is valid up to 



fourth order corrections in 1/d (see discussion in sect. 5.1). 

These results are better understood with reference to Fig.2, where the relation 
between J and J^ is plotted in the weak coupling regime (curve (a)) and in the 
strong coupling regime (curve (b)). Notice the large overlapping region around 
J^ = 0.7, where the strong and weak coupling definitions essentially coincide. 

The Douglas-Kazakov phase transition, marked by the points whose coordinates 
are in the second and third columns of Tab.IV, gives for each dimension d the 
extreme of the curve (a) below which the weak coupling solution is not anymore 
valid. For very large d this region goes as far down as J^ = 0.1 into what one would 
naively expect to be a strong coupling regime. Indeed, as already remarked, it is 
really at large d that the weak coupling solution can be already trusted just above 
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the Douglas-Kazakov transition, as the radius of the Wigner distribution is then 
small, and the approximations made in neglecting quartic terms in 9's are reliable. 

The instability of the symmetric vacuum on the other hand, marks the end of the 
region where the strong coupling expansion can be relied upon. Beyond that point 
in fact the system is forced into a configuration where presumably all the Fourier 
modes Xn of the eigenvalue distribution are excited (see discussion at the end of 
sect. 5.1) and hence non perturbative effects in J take place. For lower dimensions 
(e.g. (i = 2, 3, 4) there is an interval of values of J^ for which both curves are 
admitted. In that range the system has two classically stable configurations, one 
corresponding to the unbroken phase and one to the broken phase. 

As J^ increases the system evolves along the curve (b), then before it reaches 
the point of instability for the symmetric vacuum it jumps onto the curve (a) above 
the Douglas-Kazakov point. The exact point of this transition has been estimated 
(earlier in this section) in the fourth order strong coupling expansion and it has 
been marked in Fig. 2 for d = 2 and d = ?> with respectively a square and a triangle. 

For large values of d there is no overlap between the curves (a) and (b). This 
points to the existence of an intermediate phase where the symmetry is already 
broken, but the eigenvalue distribution is different from the semicircular one. 

6 Conclusions 

Let us now try to compare our predictions on the critical deconfinement temperature 
with the known results obtained with Montecarlo simulations. Let us stress again 
that we must expect two kinds of systematic deviations. The first one is due to the 
large N approximation. Indeed, the values that we have obtained in the large N 
limit should be considered as an upper bound, and in the cases N = 2 and A^ = 3, in 
which most of the simulations have been performed, rather large deviations should 
be expected. However, this is not an unsolvable problem. The techniques described 
in this paper can in fact be used also at finite A^, and work is in progress to do so. 
Preliminary numerical tests show that the deviations in the predicted value of the 
critical temperature are exactly in the desired direction. 

The second problem is that, in order to solve the model, we had to decouple 
space-like and time-like plaquettes in the original action. A direct consequence of 
this approximation is that the rule with which the coupling f3 changes as a function 
of Nt is essentially that of the strong coupling expansion. Once the spacelike plaque- 
ttes are decoupled, the theory becomes effectively two dimensional from the point 
of view of the strong coupling expansion, hence there is in principle no obstruction 
to extend the expansion to arbitrarly high values of f3 and At, and the continuum 
limit can be taken. The scaling law of A^ in this limit is completely fixed by eq. 
(|TT|) if we use the Wilson regularization, or equivalently by eq. (|r^) in the Heat 
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0.2 0.4 0.6 



1.2 1.4 J^ 



Fig. 2. Redefinitions of the Heat Kernel inverse coupling J^ necessary to make 
contact with the Wilson case (a) in the weak coupling regime (large J^ ), and (b) in 
the strong coupling one (small J^ )■ In the region 0.4 < J^ < 1 there is substantial 
agreement of the two redefinitions. 
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Kernel case. We see that [3/ [3 = l/Nt, or, equivalently, that 

T.= | . (81) 

with k a constant to be determined. As a direct consequence of the decouphng 
approximation, this scahng law holds in any dimensions and for any gauge group. 
The relevant point is that, in the (2+1) dimensional case, eq. ( ^l]) in fact coincides 
with the correct scaling behaviour of the whole theory. This is probably due to 
the fact that in this case the decoupled spatial theory is a two dimensional gauge 
model, and, as such, essentially trivial. Notice that a similarly good behaviour of 



the decoupling approximation in (2+1) dimensions was already noticed in [|T5[ in 
the complementary problem of describing the spacelike string tension. 

The second important point is that our value of the deconfinement temperature 
(which is effectively extracted as if we had Nt = 1) must be, already at Nt = 1, 
in the weak coupling regime of the theory, unless we would have an additional 
uncertainty in making the extrapolation to /? — > cxd. Remarkably enough, in (2 + 1) 
dimensions the values that we find fulfills also this second condition, both in the 
Wilson and in the Heat Kernel case. Notice however that these values are at the 
border of this region, a fact that will require some further caution, and will be 
better discussed in the following subsection. 

This must be contrasted with the (3+1) dimensional case. First, in the full 
theory a completely different scaling law is expected. Second, the critical point, 
at least in the Wilson case, is below the threshold of validity of the weak coupling 
expansion (see Tab. IV). These are actually two aspects of the same problem, and 
they are warning us that the result obtained within our approach must be definitely 
considered as strong coupling result. These considerations show that a more refined 
analysis, for instance by using the techniques of the (hot, twisted) Eguchi-Kawai 
model, would be important in this case. 

It should be noticed, however, that an incorrect scaling of dimensional quantities 
does not necessarily imply a total loss of predictive power. In fact, as recently 



noticed in p9[, while even the most precise and recent simulations of SU(2) and 



SU(3) models in (3+1) dimensions still show strong deviations from asymptotic 
scaling, scaling of dimensionless ratios is fulfilled even at small couplings, and very 
small values of Nt (down to Nt = 3). This opens the way to a possible improvement 
in our predictions, if other dimensional quantities besides the critical temperature, 
say the string tension, could be evaluated, for instance in the framework outlined 
in 
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Let us now discuss, in turn, the (2 + 1) and the (3 + 1) dimensional cases. 

6.1 2+1 dimensions 

In (2+1) dimensions our best prediction for the large A^ value of the critical coupling 
Jc = Pc/N'^ is Jc = 0.601. Remarkably enough, this result follows both from the 
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Wilson and the Heat Kernel regularizations, once the relation between Wilson and 
Heat Kernel couphngs is used. The result should be considered as the fourth order 
improved version of the standard mean-field result Jc = 1/d = 0.5. Comparing 
with the second and third order results we see that this should be considered as an 
upper bound, and that further orders would probably decrease the value. Keeping 
in mind the various factors of A^, we obtain 

Looking at Tab. IV and Fig. 2 we see that this value is at the border of the 
weak coupling region, hence its extrapolation toward the continuum limit {Nt — > 
oo, /3c — *■ oo) could be affected by systematic errors. A simple way to estimate these 
possible deviations is to notice that, had we chosen from the beginning a Heat 
Kernel regularization for the whole theory, then, as discussed in the introduction, 
the scahng law eq. (RTT) would have been exact. In that case we could not have used 



eq. (|61|) to map the Heat Kernel coupling to the Wilson one; instead, we would 
have had to keep the true Heat Kernel critical coupling, which turns out to be (see 
sect. 5 and Tab. IV) J^ = 0.416. Since the simulations with which we shall compare 
our results were performed with the Wilson action we shall use Jc = 0.601 in our 
comparison, but the difference between Jc and J^ gives us an idea of size of the 
possible systematic errors involved in the continuum limit. 

To the best of our knowledge, Montecarlo estimates of the deconfinement tem- 
perature in (2-1-1) dimensions only exist for the N = 2 and A^ = 3 models. 

In the case of SU(2), very precise and careful estimates of the critical temper- 
ature exist [|16| for the lattice sizes 2, 3, 4. The values are reported in Tab. V, 
together with an estimate of the continuum limit (/? -^ oo), using only the leading 
0(l//9) corrections (see [^ for details). While the deviations from scaling are still 
rather strong in the Tc/g"^ data they are, as expected, much more under control 
for the dimensionless ratio Td \fo^ whose continuum limit requires the inclusion of 
only 0(1//?^) corrections. 



A^ 


Nt 


Tc/g' 


Tj^icr) 


2 


2 


~ 0.433 


~ 1.015 


2 


3 


~ 0.417 


~ 1.063 


2 


4 


~ 0.412 


~ 1.088 


2 


oo 


0.385(10) 


1.121(8) 


3 


2 


~ 0.454 




oo [|15|| 




0.477 


0.977 


oo, Jc 




0.601 




oo, Jf 




0.416 





Tab.V. The deconfinement temperature Tc and the dimensionless ratio Td \fo as 
a junction of the lattice size in the t direction Nt in the (2+1) dimensional SU(2) 
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and SU(3) LGT. Taken from / |7^/ and jj^/ . In the sixth row we report the estimate 
of ref. ITBJ , while our predictions are reported in the last two rows. 



In the case of SU(3) there is, to the best of our knowledge, only one numerical 
result. It is the critical coupling for the Nf = 2 lattice, which turns out to be 
/3c{Nt = 2) ~ 8.17 ig. This corresponds to Tc/g^ ~ 0.454. No extrapolation to 



the continuum limit is possible in this case, however it is interesting to see that the 
value is definitely higher than those of the N = 2 case. This is a signature of the 
fact that as N increases the T^/g"^ values moves toward our large A^ limit. 

Finally, let us mention the result obtained in [jl5| (also reported in Tab. V) by 
using as additional input the dimensionless ratio a = 7iT^/3, which comes from 
the effective string approach to the interquark potential in LGT's. The resulting 
prediction for the critical temperature is 

3(iV2 - 1) 

which means, in the large N limit, T^/g"^ = 0.477, a result which again is in the 
same region of our large A^ estimates. 

6.2 3+1 dimensions 

In this case we can test our results with a much larger sample of data. The main 
problem is obviously that the scaling laws (namely the dependence of f3c on the value 
of A^j) of our simplified model and of the whole theory are completely different. This 
essentially means that, due to the decoupling of the space-like plaquettes (and to 
the fact that in this case the purely spacelike theory is non trivial), the lattice 
spacing in the time direction, with which we measure Nt in our decoupled model, is 
related in a non trivial way to the lattice spacing of the starting theory. The only 
way to use our results is then to compare our decoupled theory and the Montecarlo 
simulations of the full theory for small values of Nt, where Jc is small and we can 
trust our strong coupling expansion. In this case we must definitely study the 
Wilson action only, and no universality argument can be used. The scaling law is 
given by eq. ([T^ ) 

Our best estimate for the critical coupling in (3 + 1) dimensions is Jc = 0.378. 



Very precise data on SU(2) and SU(3) in the range A^f = 2 — 16 can be found in ||29 
Their values in terms of the A scale defined as 

together with the corresponding values of Jc = (3c/ N"^ are shown in Tables Via and 
VIb. 
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Nt 


Pc/N' 


Te/A 


2 


0.4700(8) 


29.7(2) 


3 


0.5442(8) 


41.4(3) 


4 


0.5746(2) 


42.1(1) 


5 


0.5932(11) 


40.6(5) 


6 


0.6066(8) 


38.7(3) 


8 


0.6279(10) 


36.0(4) 


16 


0.6849(25) 


32.0(8) 



Tab. Via. The critical coupling (3c/ N"^ and the corresponding deconfinement 
temperature T^/A as a function of the lattice size in the t direction Nt in the 
(3+1) dimensional SU(2) LGT. The data are taken from 



Nt 


(3c/N^ 


Te/A 


2 


~ 0.568 


-78.8 


3 


0.6167(11) 


85.7(1.0) 


4 


0.63250(2) 


75.41(2) 


6 


0.65490(6) 


63.05(4) 


8 


0.667(3) 


53.34(1.5) 


10 


0.6844(8) 


51.05(40) 


12 


0.6964(13) 


48.05(65) 


14 


0.7092(11) 


46.9(5) 


16 


0.716(6) 


44.3(2.5) 



Tab.VIb. The critical coupling fSc/N"^ and the corresponding deconfinement 
temperature Tc/A as a function of the lattice size in the t direction Nt in the 



(3+1) dimensional SU(3) LGT. The data for Nt >3 are taken from /|^/, the 
single one at Nt = 2 is taken from ^3l\] . 



The values of Jc from Tables Via and VIb corresponding to small values of 
Nt, together with some results for a large value of A^ (A^ = 81) obtained through 



Montecarlo simulations p2[ of the hot twisted Eguchi-Kawai model are compared 
with our predictions in Table VII. Our values are obtained using eq.(|T2]) and the 
value {Nt = 1, Jc = 0.378) as input. 
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Nt 


N = 2 


A^ = 3 


A^ = 81 


A^ = oo 


2 


0.4700(8) 


~ 0.568 




0.615 


3 


0.5442(8) 


0.6167(11) 




0.723 


4 


0.5746(2) 


0.63250(2) 




0.784 


5 


0.5932(11) 






0.823 


5.14 






0.695(5) 


0.828 


6 


0.6066(8) 


0.65490(6) 


0.705(5) 


0.850 



Tab. VII. The critical coupling (3c/ N"^ as a function of Nt and N . In the last 

column our predictions, according to the scaling law eq.^T^) and [3c/ N"^ = 0.378 for 

Nt = 1. The data for N = 2 and 3 are taken form Tab. Via and VIb. Those for 

N = 81 are taken from ^3^ . 



We may observe that 

a] the scahng law eq.(0) is definitely different from the scaling behaviour eq. (^4D . 
However the Montecarlo estimates of the critical couplings seem to be still far 
from the asymptotic scaling region and the values of Tc/A keep decreasing as 
Pc and Nt increase. It is interesting to notice that the critical couplings for 
small values of Nt have a dependence from Nt which is not too different from 
our strong coupling expectation eg . ([12D . 



b] As expected, the best agreement between data and predictions is obtained for 

our lowest value of A^^: Nt = 2. 

c] At fixed Nt, as N increases also Jc systematically increases, and our large A^ 

estimates seem to be the upper limit of this behaviour. 



As we have seen, our results, obtained analytically in the large A^ limit, show 
a more than qualitative agreement with Montecarlo simulations performed on the 
full theory; having used two different lattice regularizations gives us some control 
on the systematic errors associated with the continuum limit; the uncertainties due 
to the large A^ approximation and to the decoupling of the space-like degrees of 
freedom can be further reduced by studying the model at finite A^, and by using the 
Eguchi-Kawai approach. We believe that a combination of these methods may allow 
in future further improvements in the analytic determination of the deconfinement 
critical temperature. 



Appendix 



34 



We are interested in the large-N limit of the Heat Kernel on the cylinder 

K.2{gi,g2' ; J"") = Y.^'^^Xr{9i)xr{92^) (A.l) 

R 

In sec. (3) we utilized this kernel as a building block for the effective theory of 
Polyakov loops. We were therefore with the case gi = 92 = V, TrK being the 



Polyakov loop and we gave in eg. ([77| - |80D an expansion of the free energy corre- 
sponding to the partition function ( |A.1| ) in powers of 2e~ 2W , up to 4**^ order. 

We consider here the general case, sketching the procedure to obtain explicitly 
such an expansion and discussing a little its meaning. 

The kernel ( |A.1| ) is a particular case of the QCD2 partition function on manifolds 
A^G,b of genus G with b boundary components, the case A1o,2- As it is well known 



after [Q the character expansion of the QCD2 partition function on Aicb can be 



interpreted as a weighted sum over a suitable set of maps from world-sheets of all 
possible genera to M^b, that is, as a string theory. 

The "string coupling constant" of this peculiar string theory is 1/A^ - there's a 
factor of N'^~'^3 j^ front of the contributions of the maps from genus g world-sheet. 
Since we confine ourself to the N = 00 limit, i.e. we do not care about 0{j^) 
corrections in eq. ([75|), the expansion we are looking for corresponds to consider 
just the string tree level {g = 0) contributions. 

The string interpretation of the QCD2 partition function arises technically lim- 
iting the sum over the SU{N) representations to the representations obtained by 
symmetrizing tensor products of the fundamental representation and of its conju- 
gate. As shown by the Douglas-Kazakov exact solution at A^ = 00 in the case of the 
sphere this is the correct way to carry out the sum in the strong-coupling (small- 
area) regime. We too will therefore sum in eq. ( |A.1| ) over these representations. 



They can be in general labeled by means of two sets of integers {/} and {m} 
satisfying 

h>l2>...>lp>0>mi>m2>...>mg (A.2) 

The {/} and the negative of the {m} identify two Young tableaux, L and M, of 
order I = J2ih and m = — X^i'^i! the representation is then called a "composite 
representation" in [Q and denoted as (ML). Its order can be defined as / -|- m. 



When the {m} are all zero, the Casimir is given by 



di^^ = lN + l + Y.k{h-2z)-L . (A.3) 



For the representation (ML) the Casimir is 



C,'L, = Cf+C£' + 2^. (A.4) 
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Equations ( |A.3| , |A.4D are the key point allowing to obtain an expansion in powers of 



e ^/^"^ for the kernel ( |A.1| ). Indeed eq. ( [A.l|) becomes 



l-\-m 



!^2{9l,92^; J^) = Y. ^ ^^ J2 J2 X{ML){9l)X{ML){92^) 

l,m=0 LeYi MeY,^ 



2Jt 

e 



,2 



i + V, i,-(i,--2i) 



(A.5) 



Yi denotes here the set of Young tableaux made of / boxes. 

The characters X{ML){9i) can be handled in the large-N limit rewriting them in 
terms of the quantities 

Pk = j^TTg', ; a, = lTr^2' (A.6) 

and of their conjugates pl{= -^Tt {91)'^) and a^. These quantities remain finite as 
A^ —* 00, turning into the Fourier modes of the eigenvalue distributions p{6) and 
a {(f)) corresponding to gi or g2: 

Pk= fd9pie)e-'^' (A.7) 



This rewriting of the characters is accomplished by means of the Frobenius rep- 
resentation, extended to accomodate the case of composite representations (ML) 



Having made explicit all the dependences from A^ it is possible to expand explic- 
itly eq. ( [A.5|) up to a certain power (call it P) of e~2^, which means to consider 



all the representations (ML) such that / + m < P, to write the resulting expression 
as the exponential of a free energy and to neglect the subleading terms as A^ — i> 00 
in the latter, obtaining finally an expression of the form (|76[) . The expansion of 
the free energy can of course be obtained by subtracting from the expansion ( [A.5| ) 



the "disconnected" contributions (for instance at order e~7^ the contributions of 
products two rep.s of order I + m = 1). 

Expanding the last exponentials in eq. ( |A.5| ) on gets power series in 1/J^N 
and l/J^N"^. However, when extracting the 0(A^^) terms only in the free energy, 
at fixed order / + m only a polynomial in 1/J^ survive f\ The resulting expression 
makes thus sense as a strong coupling expansion. 

Now we give the free energy, eq. (|76| ) up to 3^'^ order. We use the trigonometric 
notation for the complex modes of the eigenvalue distributions p and a: 

Pk = rue'"" ; a^ = Skc'^'' ; (r^, s^ G R). (A.8) 



^This fact has a nice interpretation in the QCD2 string language; it is related to the existence 
of a limiting value for the Euler characteristic of the covering space at fixed characteristic of the 
target. 
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-1. 7^^ 



F (91,92 ;J 



e ' ' riSicos(ai -7i) + ^ — — 



-r\ - s\- -^r2s\ cos(a2 - 271 



— Tff '^1-^2 cos(2ai - 72) - -T^r\s\ cos(ai - 71) + r2S2 cos^(a2 - 72) 
+^^^^r?^?cos(2(ai-7i)) 



+ 



e 



-l/2J^\3 
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-4-6J^ , , „ . 3 



—3 rir2Si cos(ai + 02 - 37i) + -^^^.s^ cos(a3 - 871^ 



3 . 1 + J^ 
+ (6risi - —^{rlsi + ris\) jj^rlsl) cos(ai - 71) 

H -jfi rls\ cos(3(«i - 7i)) + (— ^rir2Si + — griraSi) cos(ai - a2 + 7i) 

_4 _ QjH 12 + 6 J^ 
H 7p — ?^?SiS2 cos(3ai - 71 - 72) H -^ — '^i^2SiS2 cos(ai + 02 - 7i - 72) 

n on 

-jj^r3SiS2 cos(a3 - 71 - 72) + {jjjriSiS2 + — p^?SiS2) cos(ai + 71 - 72) 

— m'^i'^2SiS2 cos(ai - a2 - 7i + 72) + —j^rls-s cos(3ai - 73) 

— Tff ^^1^253 cos(ai + a2 - 73) + 2r3S3 cos(a3 - 73) (A. 9) 

In the case 91 = 92 the free energy up to 4**^ order has been given in eq. (pO|). 
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